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Abstract: Using the stationary lightcone perturbation theory, we propose the complete 
and careful derivation the JIMWLK equation. We show that the rigorous treatment re- 
quires the knowledge of a boosted wave function with second order accuracy. Previous 
wave function approaches are incomplete and implicitly used the time ordered perturba- 
tion theory, which requires a usage of an external target field. 



The JIMWLK evolution equation is the important part of the high energy QCD. Origi- 
nally, the equation was derived using the path integral language [jl], [2], ^] (see also Q for 
generalization). The clear interpretation of the equation in terms of parton wave func- 
tion and its evolution was proposed in ||. The essential advantage of the wave function 
derivation is the explicit separation between scattering and evolution. This means that the 
projectile evolution can be treated irrespective to a target. Evolution and scattering are 
distinct and separated parts of the theory. A scattering is Wilson lines and an evolution is 
a Hamiltonian-driven task. Then, in order to calculate amplitude of a realistic process, it 
is necessary to calculate a quasiclassical gauge field generated by projectile and target. 

The main aim of this paper is to finish the incomplete derivation which was proposed 
in H where the author used a so called "ingoing" and "outgoing" wave functions. Such 
terms do not exist in the framework of the light cone perturbation theory (LCPT). Instead, 
they are elements of the old-fashion time ordered perturbation theory. Similar situation 
was observed in § where different cuts of time ordered diagrams correspond to terms of 
the JIMWLK and BK equations. Although the time ordered perturbation theory also can 
be used for JIMWLK derivation, but this requires the explicit usage of an external field, 
which plays role of time depended potential, and the special technique of time ordered 
perturbation theory, which strongly differs with technique that was used in 0. Summing 
our remarks, we can say that the work || contains improper mixing of the two perturbation 
methods which have different areas of applicability. Neither the time ordered nor stationary 
perturbation theories were used. The time ordered perturbation theory can be applied for 
the evolution of scattering amplitude of a projectile in an external field and it is hard to 
apply it to the evolution of a wave function. The wave function evolution contains more 
information, since it contains information about rapidity of partons. So the method of 
stationary perturbation theory is more general. In this paper we show that the direct usage 
of LCPT requires knowledge of the boosted wave function up to second order accuracy of 
the perturbation theory. The calculation of the second order correction is not trivial and 
requires the careful treatment with color indexes of projectile partons. The usage of pure 
LCPT allows us to present mathematically complete and rigorous derivation, which allows 
to study only projectile without any mention about a target. Currently, there is no such 
derivation in literature. We expect that our derivation is more useful for a first study of 
JIMWLK. 

We take lightcone coordinates where the coordinate x + plays role of time in the canon- 
ical quantization scheme. In order to make all constraints to second class we impose the 
following common gauge fixing conditions: 

A+ = 
d^A* = d-A~ + d i A i = 

So the gauge field A^ has only two independent degrees of freedom. After the quantization 
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(1) 



we have 

A£{$, x~)= [ -L= ( at n e *^k + x--ik3 + fi -ifc+x-+ife , j ^£ 

k+>0 



(2) 



where i = 1,2 enumerates the gluon polarizations, ef(A;) are polarization vectors which 
obey: 

^ + = 

fc + £~ - fee; = (3) 
£ T( k ) = 

We select real and orthonormal ef. 

EiEj = 8ij (4) 

Useful completeness property is 

^2 (£iu)(eiv) = uv (5) 

i=l,2 

In the one gluon Hilbert space with the basis \k, k + ) we shall use basis of transverse 
coordinates. 

(x,k + \y,p + ) = 5 2 (x- y)5{k + - p + ) (7) 
The gluon creation operator in the given transverse point is 

«U + = I e^a\^L= (8) 

Instead of momentum k + we shall often use rapidity. A momentum p + has rapidity Y 
relatively k + if 

p+ = k + e Y (9) 

The measure is 

If we increase on 5Y the rapidity of the state \P), then the wave function will be changed. 
New partons will arise in the new available rapidity window. In spirit of the RG flow we 
assume that new partons have rapidity much less than primary fast partons. 

Consider a composite dilute projectile \P). We assume that it is moving to the right. 
It scatters on a leftmoving dense target state \T). The both two considered states are 
composed from quarks and glouns. In the eikonal approximation every fast parton of \P) 
scatters independently and quasielastic S-matrix is 

+ °° 

ig J a a (x,x + )p a {x)dx + d 2 x 

S = Pe -°° (11) 



- 2 - 



where a a (x, x + ) = A^ a (x, x + ) is a classical field generated by target and p a is the color 
charge density operator. We wish to know how the scattering amplitude will be changed 
when we go to the high rapidity 

= (Psy\S\P 5 y) - (P\S\P) (12) 

where S[a] = (P\S\P). The key task is the calculation of state \P$y)- It is necessary to 
calculate this state with second order accuracy of the QCD perturbation theory. 

\P SY ) = |p(o)) + |pW) + |p(2)) (13) 

= (P«|5|P«) + (PW\S\P^) + (P^\S\PM) (14) 

The operator S does not change a parton number and parton coordinates. It acts only 
on color indexes. Hence, it is necessary to know the part of \P^) without additional soft 
gluons. 

The lightcone QCD Hamiltonian has complicated structure ||. In our calculations 
there is only one relevant interaction term. It contains a gluon emission from a color 
source. 

V = - g j A1(x,x-)J%{x,x-)d 2 xdx~ (15) 

where the current J a consists of both quarks and gluons. A slow gluon has momentum 
k + < p + , \k\ <C p + where p + is the momentum of the source parton. So only matrix 
elements like (P| J^|P) has nonzero values. The other components Ja are kinematically 
suppressed. The operator A a _ (x, x~) creates a soft gluon which has a momentum k + <C p + - 
So in (|l5| ) we can neglect x~ dependence in A°L and set x~ = 0. 

V = -g J A-(x,0)p a (x)d 2 x (16) 

p a (x) = J J+@,x-)dx- (17) 

Consider the eigenfunctions of the free Hamiltonian l^n ^)- Applying the direct analog of 
the nonrelativistic quantum mechanic perturbation theory, we find the first correction 

i^ 1) ) = ^i ) i^ 0) ) (is) 



°nfe - (0) „(0) ^ iy ^ 



,(1) _ Vfcn 
P (0) 

In our case the momentum k~ plays role of energy. Hence, energy of an emitted slow gluon 
is much larger than energy of the state |P) due to k + <ti p + and the following identity for 
massless partons 
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So the expression ( |18|) can be rewritten as 

ipW) = -if-V|p(°)) = /" !£W^L a t p"mp(°)\d 2 kdk + d 2 x (21) 

J k 2 V2k+^/{2^ M+V 

where the integration on dk + must be performed on the new opened rapidity window. 
Going to the coordinate representation, we obtain 



ipW) = / 2 g^ e <*(iT-a) 6 t . p°(£)\ P W)d 2 kdk + d 2 xd 2 y (22) 

' 7 7 ^21+7(2^)5 A;2 Rfc+.t.a'' ^ ' ^ > 

It is necessary to calculate the integral 1(f) = J kk~ 2 e lkr d 2 k. The value of the integral is 
a vector so it must has form fF(r), where F = r~ 2 (rl). 

1(f) =' ^ e ikr d 2 k = (23) 
r 2 J %i r 2 

oo 

where we used the integral f e lt dt = i . Finally, for the first order we have 

o 

|P<1,) = sJtF / '^I0^ n ^ ))dk+dxdy (24) 

This formula can be easily understood. Each projectile parton with coordinates x can 
emit a soft gluon into an arbitrary point y. The operator p a (x) measures partons in the 
projectile. The integration on x and y enumerates all possible emission possibilities. 

In the second order we want to know the state |P^ 2 )) projected on subspace without 
additional soft gluons. We can not obtain an analog of the formula (JO]) because the energy 
difference tends to zero. This lack of knowledge is compensated by the normalization 
requirement of a wave function. In our case there is an additional complication. The 
sector without soft glouns contains not only state |p(°)) but also state |p(°)) with arbitrary 
rotations of partons color indexes. So |P^ 2 )) can differ from |p(°)) on a color rotation. 

Let the |pl ) be an orthonormal basis in the space generated by arbitrary transforma- 
tions of parton color indexes in |p(°)). Let (Pn 1 '*) and |P^) be first and second corrections 

(2) 

respectively. In order to find |P„ ) we require 

+ <4Pl + <^ 2) l) (l^ 0) ) + l^) + l^ 2) >) (25) 
which in the second order gives 

(P^) + (P^) + <P{P\PP) = (26) 
Let |P1 2) ) = C nk \P^) and F mn = (P^\P^), then we have 

Cnm C mn F mn (27) 

It is clear from ( p4| ) that the matrix F mn is hermitian. So the antihermitian part of the 
matrix C nm can be arbitrary. But it does not give a contribution to the physical S-matrix. 
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The situation is the consequence of unobservability of an absolute phase of a wave function. 
We assume that the matrix C nm is hermitian. So, we have 



\p£ ) ) = -l\p£ ) )(p£W ) ) 



(28) 



The basis \Pm ) is complete so, using (24), we have the final expression for the second 
order correction to the boosted projectile wave function. 



|p(3)) 



-5Y- 



2(2tt) 3 J (z-y) 2 (z-x) 



(29) 



zxy 



where 5Y comes from the integration on the soft gluon momentum k + in the {Pm^\Pn^) 
The (P^\S\P^) term in ( |i~4|) can be evaluated with the help of the following identity. 



V ba (z)S 



(30) 



where V ba {z) is the one gluon scattering amplitude which is Wilson line in the adjoint 
representation. 



V(z) 



Pe 



ig f a a (z,z+)T% D dz+ 



(31) 



and where T\ D are the generators of the gauge group in the adjoint representation. Sub- 
stituting and |P( 2 )) into (0), we obtain the JIMWLK equation 



dS[a] 
dY 



(2vr) 3 



K Z xy {P\ 



zxy 



( -Sp a (y)p a (x) 
-p a (y)p a (x)S 
\+2V ba (z)p b (y)Sp a (x) 



\P) 



(z-y)(z- x) 



(32) 



(33) 



zxy ^_^ 2 ^__^ 2 

The equation ( |32[) can be expressed in more familiar form. We define the left and right 
derivative operators J"(x) and J°l(x) of color rotations 

1 5 



r ± { x ) 



dS_ 

dY 



(2vr) J 



K 



zxy 



ig 6a a (x, ±oo) 
-Jl(x)JUy) - J a -{x)J-(y) + 2V ba (z)J b + (x)J a _(y) 



s 



(34) 



(35) 



zxy 



It should be stressed that it is not necessary to calculate the higher orders of the lightcone 
perturbation theory with aim to improve the JIMWLK equation. All higher terms contain 
higher powers of 5Y and does not give a contribution to the differential equation. The 
corrections come only from the matrix elements of the interaction part of the lightcone 
QCD Hamiltonian (15). From the same reasons we can conclude that an emission of two or 
more gluons into the new opened phase space does not contribute to the evolution equation. 
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In the path integral language the restriction only to one soft gluon emission corresponds 
to exactness of the saddle point approximation which used for the integration of quantum 
fluctuations in the JIMWLK derivation. One trouble only with the exact evaluation of 
the saddle point. Note that all previous speculations are valid only if we study the elastic 
or quasielastic processes for which the JIMWLK has been initially intended. For inelastic 
inclusive processes we really need to calculate the higher orders perturbation theory terms. 

The lightcone QCD Hamiltonian is plagued by the existence of zero modes and requires 
a careful gauge fixing procedure. One method of studying of all orders was proposed in 
@, |§| and in § where path integrals are used. In the work ||] was considered the first 
correction to the JIMWLK equation in the dilute regime. In general we can say that the 
higher corrections gives so called pomeron loops which are very intensively studied last 
decade. 

Acknowledgment: We thank N.V. Prikhodko for feedback and useful remarks. 
References 

[1] J. Jalilian-Marian, A. Kovner, A. Leonidov and H. Weigert, Nucl. Phys. B504 (1997) 415; 
Phys. Rev. D59 (1999) 014014; J. Jalilian-Marian, A. Kovner and H. Weigert, Phys. Rev. 
D59 (1999) 014015; H. Weigert, Nucl. Phys. A 703 (2002) 823. 

[2] E.Iancu, A. Leonidov, L. McLerran, Nucl.Phys. A692 (2001) 583-645, arXiv:hep-ph/0011241. 

[3] A.H. Mueller, Phys.Lett. B523 (2001) 243-248, arXiv:hep-ph/0110169. 

[4] Y. Hatta, E. Iancu, L. McLerran, A. Stasto, D.N. Triantafyllopoulos, Nucl.Phys. A764 (2006) 
423-459, arXiv:hep-ph/0504182. 

[5] A. Kovner, arXiv:hep-ph/0508232. 

[6] S. J. Brodsky, H. C. Pauli, S. S. Pinsky, Phys. Rept. 301, 299 (1998). 

[7] A. Kovner, M. Lublinsky, Phys.Rev. D71 (2005) 085004, arXiv:hep-ph/0501198. 

[8] A. Kovner, M. Lublinsky, U. Wiedemann, arXiv:0705.1713 

[9] J.L. Albacete, N. Armesto, J.G. Milhano, JHEP 0611 (2006) 074, arXiv:hep-ph/0608095. 



-6- 



